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In this paper, we investigate tlie atom-molecule conversion dynamics of a generalized many-body 
model that includes the atom-atom, atom-molecule, and molecule-molecule interactions, empha- 
sizing the efficiency of the Feshbach molecular formation. We show that the picture of two-body 
molecular production depicted by the Landau-Zener model is significantly altered: The energy levels 
are dramatically distorted and the conversion efficiency is suppressed by the particle interactions. 
According to the rule of constant action and with the help of phase-space analysis, we derive an 
analytical expression for the conversion efficiency in the adiabatic limit. It shows a ceiling for the 
conversion efficiency when the interaction strength is larger than a critical value. We further derive 
a closed equation for the conversion efficiency with the stationary phase approximation. In the sud- 
den limit, the conversion efficiency is twice that predicted by the two-body Landau-Zener formula. 
Our analytical formula has been confirmed by numerical calculations. 
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I. INTRODUCTION 

The production of ultracold diatomic molecules in 
bosonic systems is an exciting area of research with 
important appUcations ranging from the production of 
molecular Bose-Einstein condensates (BECs)[l] to the 
study of chemical reaction dynamics [3] ■ A widely used 
production technique involves the association of ultracold 
atoms into very weakly bound diatomic molecules by ap- 
plying a time varying magnetic field in the vicinity of a 
Feshbach resonance[3l|J|. The underlying conversion dy- 
namics are usually described by the Landau-Zener (LZ) 
model[^]. In this model, the Feshbach molecular produc- 
tion is discussed under a two-body configuration where 
a single pair of atoms is converted to a molecule at an 
avoided-crossing between atomic energy level and molec- 
ular energy level while the molecular energy is lifted by 
an applied linearly sweeping magnetic field. Thus, the 
molecular production efficiency is expected to be an ex- 
ponential Landau-Zener typejy, |7|. 

However, in the above two-body modelsjy, 0] and 
their many-body extension[8|, |9|, interactions between 
particles such as the atom-atom, atom-molecule, and 
molecule-molecule interactions, were not taken into ac- 
count. Recently, Santos et. al. investigated a general 
many-body two channel model in which the above par- 
ticle interactions have been taken into account in the 
Hamiltonian [10| . More interestingly, the inclusion of 
these particle interactions brings richer dynamics to the 
Josephson oscillation of the coupled atom-molecule sys- 
tem. 

In this work, we extend to study the conversion dy- 
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namics of this general many-body model to determine 
how these interactions affect molecular conversion as the 
magnetic field is swept across the Feshbach resonance. 

In the adiabatic limit, we cast the many-body model 
into an effective classical Hamiltonian in the meanfield 
approximation. Such classical analog facilitates us to an- 
alyze the conversion dynamics within the frame of classi- 
cal adiabatic theory [ll|, |l2|. In this framework the con- 
version dynamics are determined by two scaled parame- 
ters U and Q, where the nonlinear interaction parameter 
U is proportional to the particle interactions while the 
coupling parameter f2 has a close relation to resonance 
width. We find that the conversion efficiency can be up 
to 100% when the scaled nonlinear parameter is smaller 
than a critical value, i.e., U/V, < •\/2/4, while there is 
a ceiling of less than 100% on the conversion efficiency 
when U/fl > \/2/4. The ceiling of the conversion ef- 
ficiency can be calculated analytically according to the 
principle of constant action. Based on this expression, 
we can also obtain the scaling laws at the critical point. 

In the sudden limit, we derived a closed equation for 
the conversion efficiency with the stationary phase ap- 
proximation. We find that the conversion efficiency is 
twice of that predicted by a standard two-body Landau- 
Zener formula when the sweeping rate is very large. This 
many body effect is consistent with the theoretical anal- 
ysis based on the renormalized Landau-Zener formula in 
Ref. i9i1 t hat has been used to explain the experimental 
data]il|. 

Our paper is organized as follows: In Sec. II we intro- 
duce our model. In Sec. Ill we make a comprehensive 
analysis of the dynamics of Feshbach molecular forma- 
tion and show how the particle interactions affect the 
conversion efficiency. Sec. IV is our conclusion. 



II. MODEL 

A widely used molecular production technique involves 
the association of ultracold atom pairs (open channel) 
into very weakly bound diatomic molecules (close chan- 
nel) by applying a time varying magnetic field in the 
vicinity of a Feshbach resonance [3|, |J]. Precisely speak- 
ing, atoms are converted into molecules at an avoided- 
crossing between the atomic energy level and the molec- 
ular energy level while the molecular energy is lifted by 
an applied linearly sweeping magnetic field. The under- 
lying conversion dynamics are properly described by a 
two channel model. The two-channel model Hamilto- 
nian that includes the atom-atom, atom-molecule, and 
molecule- molecule interactions, takes following form [lO|: 



H = ^a^a^aa+^b^b^bb- 
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where aJ is the creation operator for an atomic mode 
while 6^ creates a molecular mode. Ca, £& are the chem- 
ical potential of the atomic mode and molecular mode, 
respectively. In experiments, the external magnetic field 
is linearly swept B{t) = Bt, and crosses the Fesh- 
bach resonance at Bq, thus 2ea — £& = fJ.co[B{t) — Bq]. 
Here, iico is the difference between the magnetic mo- 
ments of a molecule and a pair of separated atoms. 
Lu = y^AnhPabgABfico/rn denotes the amplitude for the 
interconversion of atoms and molecules due to the Fes- 
hbach resonance, in which m is the mass of a bosonic 
atom, ahg is the background scattering length, and AB 
is the width of the resonance. In contrast to the early 
Hamiltonians[^, |9|, we have included the background 
scattering interactions between atom-atom, molecule- 
molecule, and atom- molecule, denoted by Ua,Ub, and Uab, 
respectively, ui — Anfi^ai/mi {i — a,b,ab) with ai and 
rrii denoting the background scattering length between 
corresponding particles and their (reduced) mass, respec- 
tively. The whole number of the system N — Na + 2Nb is 
a constant and is very large in the present experiments. 
Here, A''^ = a^a and A'^;, = b^b are the atom number and 
molecule number, respectively. We introduce the param- 
eter V to denote the quantized volume of the trapped 
particles, therefore n = N/V is the mean density of the 
initial bosonic atoms. 

Choosing the Fock states as the basis, the Schrodinger 
equation is written as 



where 






(2) 



) = J:tloC^\N - 2z,z), |A^ - 2*,z) = 
^^^={a^a^f/'-\brr\0,Q) {^ = 0,...,7V/2) are 

Fock states, and Ci is the probability amplitudes of the 
corresponding Fock states. The normalization condi- 
tion is that ^j |cip = 1. The Hamitonian matrix el- 
ements are Hij = {N — 2i,i\H\N — 2j,j) . For i = j, 



2i)ea + isb] For i ^ j, F^ = except -ffi,i+i = i?i+i,i = 
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For the simplest case of A^ = 2, the above Schrodinger 
equation reduces to the following two-level system of 
standard Landau-Zener type. 



Co 
dt \ ci 



2A 2n 

2n -2A 



(3) 



where |cop and |cip denote the population of atoms 
and molecules, respectively. The energy bias is A = 
(2ea — Cfc + 2nua) /4 and the coupling strength is given 
by 57 = uj^/n/2. Initially, all particles populate in the 
lower level of the atomic state, i.e., cq = 1, ci = 0. When 
the external magnetic field is linearly swept across the 
Feshbach resonance at A ~ 0, a fraction of atoms will 
be converted to molecules at the avoided-crossing of en- 
ergy levels. The conversion efficiency as a function of 
the sweeping rate (i.e., a = A = ^coB/A) and coupling 
strength, takes the formQ, 



X = 1 - Tj^ = 1 - cxp(^ 
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It is interesting to notice that in the two-body molec- 
ular production picture the particle interactions do not 
affect the conversion efficiency. This is quite different 
from the many-body picture as we discuss later. 



III. MANY-BODY EFFECTS (iV -^ oo) 

As the total particle number A^ increases, Bq.© is 
no longer analytically solvable. In addition, the com- 
putational demand increases dramatically as A^ becomes 
very large. In the mean-field limit where A^ — > oo, the 
quantum fluctuation is negligible. It is appropriate to 
replace all the quantum operators with c- numbers, thus 
the Heisenberg equations are casted into the following 
nonlinear Schrodinger equation. 



dt 



= H 
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-(2ea - eb + 2nwQ - -nwb), (8) 
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and the total population is normalized to a unit, i.e., 
|a|V2|6|' = l. 
We introduce the canonical transformation, 



S 



26'a — Oil, 



(10) 

(11) 



where 9a = arga is the phase of the atomic mode and 
Oh = argfe is the phase of the molecular mode. The 
quantum system is equivalent to the following classical 
Hamiltonian, 

7^ = -2[/S'2 + 2A5 + 2f7(l + 5)\/]~^cos6'. (12) 

The new canonical variables satisfy 

dS dn 



dt 
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(13) 
(14) 



Before proceeding, we would like to draw attention to 
some interesting characteristics of our present Hamilto- 
nian. Compared to the well known nonlinear Hamil- 
tonian that describes the tunneling dynamics of BEC 
atoms in a double-well potential or between two inter- 
nal states [IJ], the nonlinearity {i.e., the Hamiltonian de- 
pends on the instantaneous wavefunction as well as its 
conjugate) of the present Hamiltonian ([6]) is much more 
complex. The nonlinearity stems both from the diago- 
nal and off-diagonal terms and is preserved even without 
taking into account the particle interactions. In addi- 
tion, the absence of hermicity as well as the lack of C/(l) 
invariance of the Hamiltonian restrict the mean-field mo- 
tion to a "tear-drop" shaped equal-single-pair-entropy 
surface I la. I Iq. rather than the surface of a Bloch sphere. 
As we will show latter, the interplay of these new features 
leads to a very different energy level structure and con- 
version dynamics. 

In the following calculations, the coupling constant CI 
is chosen as the energy scale, thus weak nonlinearity and 
strong nonlinearity refer to U/il << 1 and U/ft >> 1, 
respectively. 



A. General properties 

We first show how the nonlinear interactions lead to 
the deformation of the eigen-energy levels. The eigen- 
statcs of the system satisfy that 



H 



= /" 



a 
26 



(15) 



Solving the above nonlinear equations together with total 
particle conservation condition |ap-|-2|6p = 1, we readily 
obtain the chemical potential /i and the eigen-state (a, b). 
The eigen-energies can be derived according to the rela- 
tionship £ = n/2 + ^|6|2 + A|a|V2 -I- 4U\b\^ - 2U\a\^\b\^. 
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FIG. 1: Adiabatic energy levels for different nonlinear inter- 
action strength: (a) f/ = 0, (b)[/ = 0.2, {c)U = 2, {d)U = 30. 
In all cases, = 1, the solid lines represent stable eigen-states, 
and the dotted lines between Ai = U — \p2 and A2 = \J -\- \/2 
correspond to unstable states. When U > \/2/4, a loop struc- 
ture appears at the lower energy level. The loop expands as 
U increases. 



Their dependence on the parameters is plotted in Fig[T] 
for the cases of linearity, weak nonlinearity, and strong 
nonlinearity, respectively. In the linear case [U = 0, 
FiglIJa)], the energy levels are center symmetric. There 
are only two eigen-states when |A| is large enough, one 
for atomic mode and the other for molecular mode. 
When |A| /rj < ^/2 , there is an additional eigen-state 
with S = —1, represented by the dotted line in FiglTJ 
This eigen-state is dynamically unstable. With the ap- 
pearance of nonlinear interaction, the symmetry of the 
energy levels breaks down. For the weak nonlinear case 
U/il < \/2/4 [see FiglUb)], the energy level structure is 
very similar to that of the linear case except for a slight 
shift. However, when U/il > a/2/4, a loop structure ap- 
pears at the lower energy level. The loop expands as 
U increases, and the gap between the upper and lower 
energy level becomes narrower and narrower. Such de- 
formation of energy levels consequently leads to very dif- 
ferent conversion dynamics. 

Consider the adiabatic evolution of the system start- 
ing from the atomic mode at the left side of the lower 
energy level. When U is small, e.g., in Fig[lja), the evo- 
lution of the system follows the solid line, converting all 
atoms into molecules. However, when J7/il > \/2/4 as 
in Figlljc), the system moves steadily from the left side 
to the critical point C. After that, there is no way to go 
further except to jump to the upper and lower levels. As 
that fraction of atoms tunnels to the upper level, they 
are not converted into molecules. The situation becomes 
even worse when U is very large: the critical point is 
much closer to the upper level and far away from the 
lower one, thus the system will jump to the upper level 
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FIG. 2: (Color online) Conversion efficiency as a function of 
the sweeping rate ajQ^ for various interaction. 



more easily, see Figjljd). As a result, almost all atoms 
can not be converted into molecules. 

The above simple analysis is confirmed by our numer- 
ical results, which are plotted in Fig[21 In our calcula- 
tions, the 4-5th Runge-Kutta step-adaptive algorithm is 
used to solve the nonlinear Schrodinger equation ([5]) . The 
conversion eSiciency as a function of the sweeping rate a 
is plotted against the nonlinear parameters ranging from 
weak nonlinearity to strong nonlinearity. Fig [2] shows 
that (1) the nonlinear interaction always suppresses the 
conversion from atoms to molecules. For example, in the 
case of strong nonlinearity CZ/il — 30, the conversion ef- 
ficiency is only about 10%. (2) The conversion efficiency 
decreases monotonously as the sweeping rate increases. 
In the following sections, we will further derive analyt- 
ical formulas for the conversation efficiency in the two 
limit cases, namely, the adiabatic and sudden limit, cor- 
responding to ajVi} << 1 and ajVi} >> 1. 



B. Adiabatic limit 

In the adiabatic limit where the external field varies 
slowly compared with the intrinsic motion of the system, 
the conversion dynamics are entirely determined by the 
phase-space structure evolution of the classical Hamilto- 
nian (J12p . The fixed points (i.e., the energy extrema of 
the classical Hamiltonian) on the phase space correspond 
to the quantum eigenstates. For example, in Fig[3{d), 
the elliptical point (red circle) corresponds to the atomic 
mode at the left side of the lower energy level in Figlljc), 
and the saddle point (blue circle) corresponds to the un- 
stable eigenstate at the tip of the loop. According to 
the classical adiabatic theory, when the energy bias A 
is modulated adiabatically, a closed orbit in the phase- 




FIG. 3: (Color online) The phase space evolution as A 
changes adiabatically. The left column and right column rep- 
resent the linear and nonlinear case, respectively. The arrows 
indicate the shifting direction of the fixed points as A in- 
creases (see text for detail). 



space remains closed and the action 

I=— I Sde 
2ttJ 



(16) 



stays invariant in time. The action equals the phase- 
space area enclosed by the closed orbit, and is zero when 
the orbit shrinks to a fixed point. 

For the case of U /Vt < \/2/4, the initial state [atomic 
mode at the left side of the lower energy level in Fig[TI^b)] 
is prepared at an elliptical point on the phase space [red 
circle in FiglSja)]. The elliptical point evolves contin- 
uously from the boundary line of 5'= 1 to 5* = — 1 
throughout the entire sweeping of A. According to the 
rule of constant action, it is thus expected that the 
system follows the elliptical point and finally reaches 
5 = — 1, implying the entire conversion of atoms into 
molecules, i.e., the conversion efficiency is x = 1. 

However, for the case of [//57 > \/2/4, the elliptical 
point [red circle in Fig[3jd)] will collide with a saddle 
point [blue circle in Figl2{d)] when A = Ac [see Figl3Je)]. 
After this collision, the system enters an new obit [ma- 
genta line in Figl^l^e)] with Ti. = He, and evolves adiabati- 
cally for A > Ac according to the rule of constant action, 
which is now nonzero. This obit eventually evolves into 
a straight line of constant 5* [magenta line in FiglJI^f)]. 
With these considerations, we can obtain the conversion 
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FIG. 4: (Color online) Comparison between analytical cal- 
culation (red lines) and numerical simulation (full circles) 
in the adiabatic limit, excellent agreement is clearly seen. 

Uo/n = %/2/4. 



efficiency in the adiabatic limit, 



(17) 



To work out the explicit expression of x, we first need 
to determine the critical point C. For this purpose, we 
notice that point C (with 6 = tt) is a. double root of 
Eq.dHlllTll), thus 



89 
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2(1 - Sc^/^ 



= 0. 



(18) 



The critical energy bias Ac and obit energy Tic can be 
obtained through Eq. (fT4|) and Eg. (fT2|) . respectively 
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-2USc^ + 2AcSc - 217(1 + Sc)\/^-Sc- (20) 



Once these critical values are given, the whole orbit [ma- 
genta line in Fig[3]Je)] passing through the critical point 
can be determined by 



Tic 



2US^ - 2AcS 



coso = 



2n{i + s)VT^^ 



We calculate the action by changing the parameter of 
the integration. 



= — (fodS 
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1 r^- 



He + 2C/S'2 - 2AcS' 
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Here, 5min can be determined by Eq. ([2T|) with 9 = 0, i.e., 

He + 2U b,-,.:,. — zAcO.nin 



1 



2r!(l + S^in)VT^^Sn 



(23) 



The results based on the above analytical calculation 
are compared with those obtained from numerical simula- 
tion in FiglHand excellent agreement is obtained. When 
U/il < V2/4 the conversion efficiency is always a con- 
stant, while when U/fl > \/2/4 the conversion efficiency 
decreases suddenly. This implies a phase transition at 
the critical point of Uq/Q — \/2/4. To confirm the oc- 
currence of phase transition and determine its order, we 
first need to obtain the scaling law around the critical 
point. For this purpose, we introduce a small variable 
S = U /Q, — '\/2/4, and calculate the critical values with 
perturbation theory, 
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The critical orbit is determined by (let S = — 1 + a;) 
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^S, a;min ~ ^^^6^. The integration of 
Eq. ip^ results in a power law, 
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(28) 



Clearly, both x ^-nd its first-order derivative are con- 
tinuous at the critical point. However, its second-order 
derivative turns to be discontinuous, indicating the ap- 
pearance of a second order phase transition at the critical 
point. 

The power law at the asymptotic regime of U/fl —f cxd 
can also be obtained with the above treatment. In this 
case, we introduce a new small variable A = (U/fl)~^ 
and use this small variable to expand the critical val- 
ues, Se = l- i(2A)2/3, 5„,in = 1 - |(2A)2/3, Ajn = 
2A-1- 3|1a-i/3 + 3|Iai/3, Hc/n = 2A-i-3^A-V3 + 



(21) ^^A^/^. After integrating Eq.([22]), we finally come to a 



power law of the conversion efficiency. 
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We would like to mention that, in the nonlinear 
Landau-Zener model describing the tunneling dynamics 
of BEC atoms in a double-well potential or between two 
internal states, the power laws at the critical and asymp- 
totic regime are very different |14| . At the critical regime, 
both the prefactor and exponent of the power law is dif- 
ferent from that in Eg. ([28]) . While at the asymptotic 
regime, only the prefactor is different compared with 
Eg. I 
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FIG. 5: (Color online) Comparison between analytical calcu- 
lation (solid lines) and numerical simulation (scatters) in the 
sudden limit. Good agreement is obtained, especially for the 



rapidly sweeping case a/0. 



100. 



C. Sudden limit 

The sudden limit corresponds to nonadiabatic conver- 
sion. The conversion efficiency is not strongly related 
to the structure of the energy levels. In this limit, we 
can derive the analytical expression of the conversion effi- 
ciency using the stationary phase approximation (SPA). 
Because of the large sweeping rate a, a quantum state 
would stay on the initial level most of the time. Thus 
the amplitude b in the Schrodinger equation ([5]) remains 
small and \a\ ~ 1 all the time. A perturbation treatment 
of the problem becomes adequate. 

We begin with the variable transformation, 



A + 2C/(2|6|^-|a|^) dt\ ,{30) 



a exp < —I 



b = b expU 2A + 4C/(2 |6r - \a\) dt ^ . (31) 

As a result, the diagonal terms in the Hamiltonian are 
transformed away, and the evolution equation of b be- 
comes: 



t ^ 



^ = -2in{a)^ exp \-i I 4A + 8U{2 \b\^ - |a|'; 
at [ Jq I 



here a ^ 1 all the time and \a\ + 2\b\ = 1, thus 



b=-2inl dtcxpl-i I 4A + 8t/(4|6|^-i; 



dt 

(32) 



dt 



(33) 



We need to calculate the above integral self- 
consistently. Due to the large a, the nonlinear term in 



the exponent generally gives a rapid phase oscillation, 
which makes the integral small. The dominant contri- 
bution comes from the stationary point to of the phase 
around which we have 



AA + 8UiA\b\^ -l)^ait-to), 



with 
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dt 
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Since |6|^ 



16 P, then we have 



\b\' = i2nf 



(it exp 



-a{t-taf 



(36) 



This expression can be differentiated and evaluated at 
time to- A standard Fresnel integral with the result 
[d\b\ /dt]tg — (217) \/Tr/a is obtained. Combining this 



with relation (|35p . we come to a closed equation for a, 



a = Aa + 32U(2nf J^. 
V a 



The conversion efficiency 
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Substituting it into Eq. (|37p , we obtain the analytical 
expression for the conversation efficiency. 
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In FigO we have calculated the conversion efficiencies 
using Eg. pop and compared them with the numerical 
results obtained by directly integrating the Schrodinger 
equation ([5]), where a good agreement is shown, espe- 
cially for the rapidly sweeping case, e.g., a/Vt^ = 100. 
In this case, the second term on the right hand side of 
Eg. ([39)1 can be neglected, thus we get x = 47rJ7^/Q!. This 
result from our many-body mean-field theory is twice 
that obtained from the standard two-body Landau-Zener 
formula, i.e., the first order expansion of Eq.(|3]). This 
many body effect is consistent with the theoretical anal- 
ysis based on a renormalized Landau-Zener formula in 
Ref.[|. 



IV. CONCLUSION AND DISCUSSION 

In conclusion, we have both numerically and analyt- 
ically investigated a generalized many-body model that 
includes the atom-atom, atom-molecule, and molecule- 
molecule background scattering interactions, emphasiz- 
ing the dynamics of Feshbach molecular formation. Com- 
pared to the simple two-body molecular production pic- 
ture depicted by the Landau-Zener model, the conver- 
sion dynamics become more complicated because of the 



interplay of the many-body effect and the nonhnear in- 
teraction effect. 

In the adiabatic hmit, the nonhnear particle interac- 
tions dominate while in the sudden limit the many-body 
effect becomes significant. In the adiabatic limit, the en- 
ergy level structure is dramatically distorted due to the 
particle interactions. Especially when the nonlinear pa- 
rameter goes beyond a critical value, a loop appears at 
the tip of the lower energy level. This loop structure 
consequently leads to the break down of adiabaticicy, 
and correspondingly a ceiling of less than 100% on the 
conversion efficiency. The above striking phenomenon 
emerges when the effective interaction parameter [7/17 is 
large. Because the above parameter is proportional to the 
square root of the density and inversely proportional to 
the square root of the resonance width, the above phe- 
nomena are expected to be observable in experiments 
performed with atom clouds of relative high density and 
narrow Feshbach resonance width. 



In the sudden limit, we derived a closed equation for 
the conversion efficiency with the stationary phase ap- 
proximation. We found that the conversion efficiency is 
twice that predicted by the standard two-body Landau- 
Zener formula when the sweeping rate is very large. This 
many-body effect is consistent with the theoretical anal- 
ysis based on a renormalized Landau- Zener formula in 
Ref. [9|, that has been exploited to explain the experi- 
mental data of JILA[13] and agrees well. 
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